Observation of Electronic Ferroelectric Polarization in Multiferroic YMn205 
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We report the observation of a magnetic polarization of the O 2p-states in YMn205 through 
the use of soft X-ray resonant scattering at the oxygen K-edge. Remarkably, we find that the 
temperature dependence of the integrated intensity of this signal closely follows the macroscopic 
electric polarization, and hence is proportional to the ferroelectric order parameter. This is in 
contrast to the temperature dependence observed at the MnLs-edge, which reflects the Mn magnetic 
order parameter. First principle calculations provide a microscopic understanding of these results 
and show that a spin-dependent hybridization of O 2p- and Mn 3c/-states results in a purely electronic 
contribution to the ferroelectric polarization, which can exist in the absence of lattice distortions. 

PACS numbers: 75.85. +t, 77.80.-e, 64.70.Rh, 61.05.C- 



Multiferroic materials with a coupled ordering of elec- 
tric and magnetic moments are extremely interesting for 
technological applications. They could, for example, be 
used to build a magnetic computer memory that could 
be switched with an electrical field (see e. g. [Dill)- This 
has motivated a long history of research on these materi- 
als. However, compounds with the desired properties are 
very rare. Magnetism and ferroelectricity seldom coexist 
and even when they do, the coupling between them is 
typically extremely weak. The discovery of a number of 
transition metal oxides (TMOs) with a strong coupling 
[3, 4 has therefore generated much excitement. 

Attention has quickly turned to understanding the 
mechanism for this coupling. Here, models based on 
structural distortions, which optimize the energy of the 
magnetic system, play a central role [1-4 . Since these 
scenarios focus on magnetically driven distortions of an 
ionic lattice, we will refer to the corresponding ionic fer- 
roelectric (FE) polarization as Vion- 

Transition metal oxides are, however, notorious for 
having strong electronic correlations and interactions 
that not only involve lattice and spins, but also charges 
and orbitals. Further, it is well known that covalency 
between the transition metal ion and the oxygen plays a 
major role in most TMOs [5 . This has led to the sug- 
gestion that there may be another contribution to the 
polarization, i.e., one due to the valence electrons. Re- 
cent theoretical studies point to such an electronic con- 
tribution, Vei to the FE polarization [6l|7]- In contrast 
to T^ion, Vei is not primarily related to structural distor- 
tions, but rather to the covalency between the oxygen 
and the transition metal, i.e., it is a polarization that 
can exist even if there is no movement of the ion cores 
from their undistorted positions. Whether or not such 
a Vei is in fact significant in these materials is currently 
hotly debated. While there is some indirect evidence for 
it to date, direct experimental tests of this picture 



have been lacking. 

In this letter, we report the observation of Ve\, in mul- 
tiferroic YMn2 05 by means of resonant soft X-ray scat- 
tering (RSXS) at the O i^-edge. The RSXS signal at the 
O i^-edge is due to a redistribution of valence charge, re- 
sulting in a magnetic polarization of the O 2p-states. The 
corresponding magnetic signal at the O K-edge does not 
follow the magnetic scattering from the Mn moments, but 
rather closely follows the ferroelectric polarization. Em- 
ploying density functional theory (DFT) calculations, we 
provide a microscopic explanation of this observation and 
show that strongly spin-dependent hybridizations result 
in a purely electronic contribution Vei to the macroscopic 
FE polarization V. Our RSXS measurements at the O K- 
edge are measurements of this Vei- We emphasize that 
Vei is possible even in the absence of any lattice distor- 
tion. 
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FIG. 1. In YMn205 Mnl-octahedral chains (blue) in the c- 
direction are linked by Mn2-pyramidal dimers (brown) in the 
a6-plane. The grey spheres are Y, the red O ^ilOj. The green 
arrows represent the commensurate spin structure at 25 K as 
determined by neutron diffraction (m I12j . The dashed line 
in the center indicates the plane shown in Fig.[4] 
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FIG. 2. (a) and (b) show the resonant X-ray scattering in the 
(ifOL)-plane as measured at the Mn La-edge (644.6 eV) and 
O K-edge (529.6 eV), respectively. The observed peak moves 
from the commensurate position (1/2, 0, 1/4) (25 K, right) to 
(0.48,0,0.29) (12 K, left), in the incommensurate phase. 



The crystal structure of YMn205 is illustrated in 
Fig.[l] Manganese occupies two different crystallographic 
sites; the first coordinated by an oxygen octahedron 
(Mnl) and the second by an oxygen square-based pyra- 
mid (Mn2) [9 ,T0^. The Mnl-octahedra form edge-sharing 
chains along the c-direction that are linked within the 
a6-plane by Mn2-pyramid dimers. Within an ionic de- 
scription, the very different local coordination of Mnl 
and Mn2 results in different formal valencies of these two 
sites, namely 4+ {3d^) for Mnl and 3+ (3d^) for Mn2 
[13, 

As a function of temperature, YMn2 05 undergoes a 
sequence of phase transitions, which demonstrate the 
strong coupling between the magnetic and FE order: 
Upon cooling, a transition from paramagnetic to in- 
commensurate antiferromagnetic (HTIC-phase) occurs 
at Ti = 45K [T4 . At T2 = 39K the antiferromagnetic 
order becomes commensurate (C-phase) [14 and a finite 
V appears [15]. Upon further cooling, low-temperature 
incommensurate magnetic order sets in at T3 = 19 K 
(LTIC-phase) [14 , which is connected to a sign change 
and a reduction in magnitude of V [15] . 

The YMn2 05 single crystals used for our study were 
grown using the PbO-PbF2 flux method [16] and charac- 
terized by specific heat, magnetization and X-ray diffrac- 
tion measurements, which were all in excellent agree- 
ment with previously published results [9l [TOl [17] . For 



the RSXS experiments, a YMn205 sample with a (201)- 
surface normal was polished using 0.1 /im diamond films. 
Experiments were performed at the UE46-PGM1 beam- 
line of BESSY II at the Helmholtz-Zentrum Berlin and 
the XI A2 beamline at the NSLS, Brookhaven National 
Laboratory. The former setup realizes a horizontal scat- 
tering experiment with variable incoming polarization 
and is equipped with a photodiode point detector. The 
scattering plane of the NSLS-chamber was vertical and 
the incoming polarization horizontal. In this case, two- 
dimensional cuts in reciprocal space were recorded using 
a Princeton Instruments PI-MTE CCD. In each case, the 
sample was cooled using a liquid helium cryostat. 

In Fig.[2) RSXS-intensity distributions in the {HOL)- 
plane are presented for the LTIC- (left column) and the 
C-phase (right column). The two (i^OL)-maps in the up- 
per row of the figure were taken with the photon energy 
tuned to the main absorption peak of the Mn I/3-edge at 
644.6 eV [18]. At this photon energy RSXS becomes ex- 
tremely sensitive to magnetic order on the Mn-sites [19]. 
In the C-phase, we observed a wave vector of (1/2, 0, 1/4) 
which was in agreement with the magnetic wave-vector 
determined by neutron diffraction [14 . In the LTIC 
phase, an incommensurate wavevector (0.48, 0, 0.29) was 
observed, again in very good agreement with previous 
neutron diffraction measurements [14, 20 . 

More surprisingly, and more importantly, we also ob- 
served a superlattice reflection at the magnetic wave- 
vector when the photon energy was tuned to the O K- 
edge at 529.6 eV [21]. This is shown in Fig.[2](b). The 
transition from the C-phase into the LTIC-phase is 
also clearly seen at the O i^-edge, as demonstrated in 
Fig.[2](b). We attribute the scattering at the Oi^-edge 
to a magnetic polarization of the O 2p states. This is 
for two reasons. First, the wavevectors measured at the 
Oi^T-edge and the MnL2,3-edge are identical. Second, a 
magnetic polarization of the O 2p states was also found in 
previous measurements on the iso-structural compound 
TbMn205 [21 . X-ray magnetic circular dichroism mea- 
surements (XMCD) on related transition metal oxides 
have shown that the magnetic signal at the oxygen K- 
edge can be interpreted as an orbital moment transfered 
to the oxygen via strong hybridization between the Mn 3d 
and 2p states [22l[23], consistent with our results. 

The width of the superlattice peak along q is about 
three times larger at the Mn La-edge than at the O K- 
edge. This is simply due to the different X-ray penetra- 
tion depths at the two edges. An estimate from our XAS 
spectra and the tabulated values of Henke [24], yields 
X-ray penetration depths of 0.04 //m and 0.12 jam at the 
MnLs- and OiiT-edges, respectively, entirely consistent 
with the different observed peak widths. 

We now turn to the temperature dependences of the 
superlattice reflections as measured at the Mn La-edge 
and the Oi^-edge. These are shown in Fig.[3| In each 
case, the intensity of the reflections was accurately and 
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FIG. 3. Temperature dependence of the soft X-ray integrated 
intensity as measured at the MnLa- (a) and O i^- edges (b). 
Good agreement is seen with the reported (sohd hues) inte- 
grated intensity observed by neutron diffraction data (a) [14] 
and magnitude of the electric polarization \V\ (b) [15], re- 
spectively. The phases are low temperature incommensurate 
(LTIC) weak ferroelectric (wFE), commensurate (C) ferro- 
electric (FE) and high temperature incommensurate (HTIC) 
paraelectric (PE). 



efficiently integrated over a full 3D volume in recipro- 
cal space using the CCD camera. The integrated in- 
tensity measured at the Mn La-edge is compared to the 
corresponding magnetic neutron scattering data |^ in 
Fig.[3](a). The very good agreement between the two 
confirms that the RSXS signal at the Mn Ls-edge is mea- 
suring the square of the magnetic order parameter - just 
as neutron scattering does. 

As can be seen in Fig.[3](b), the temperature depen- 
dence observed at the Oi^T-edge is quite different and 
does not follow the magnetic order parameter. Instead, 
it drops dramatically at the C-LTIC transition and does 
not recover at low temperatures, reaching only about 20% 
of its maximum value in the commensurate phase. This 
is the central result of this study: the RSXS signal at the 
O K-edge closely follows the absolute value of the macro- 
scopic ferroelectric polarization |7^|, which is plotted in 
Fig.[3](b) [15 . Within the error of our experiment, the 
two quantities have identical temperature dependences, 
even through the transition from the C-phase into the 
LTIC-phase. This establishes a direct relation between 
the electronic modulation of the O 2p states and the fer- 
roelectricity of YMn205. 

Motivated by these results, we have performed 
LSDA+U calculations in order to clarify this relation- 



ship. The ah initio calculations were performed employ- 
ing the full potential WIEN2K code [25j, using 250 k-points 
in the irreducible Brillouin zone. We set U = 5eV and 
J = 0.88 eV (cf. Ref. 7) and performed calculations with- 
out spin-orbit interactions. The calculations were based 
on the atomic positions, which we determined by means 
of standard single crystal X-ray diffraction for our sam- 
ples in the C-phase. The resulting structural refinements 
are in very good agreement with previous studies [9l [TO] 
and gave Pbam as the space group of YMn205. As in 
previous studies, no additional structural distortion or 
symmetry reduction could be resolved in the C-phase 
[26] , Therefore we chose the maximal non-isomorphic 
subgroup of Pbam that allows for the experimentally ob- 
served V along the 6-direction, namely Pb2im [27] . 

Since the magnetic unit cell of the real material is 
large and we are only interested in qualitative effects at 
this point, we studied a simplified spin structure, which 
approximates the observed C-phase and can still be de- 
scribed using the Pb2im unit cell. In this symmetry, all 
the Mnl-sites are equivalent and were set to have parallel 
spins. In addition, there are two inequivalent Mn2-sites, 
Mn2a and Mn2b, in Pb2im. We set the spin of Mn2a 
to be parallel and that of Mn2b to be antiparallel to the 
Mnl-spin. 
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FIG. 4. Results from LSDA+U calculations in Pb2im sym- 
metry, (a) Charge density pc = p(t) + p(i) of the valence 
states within 1.5 eV of the Fermi level in the 6^c-plane. The 
6^-axis is tilted by 23.5° from the 6-axis (Fig.[l]). (b) Projec- 
tions of pc shown in (a) (blue line) and the corresponding spin 
density ps — p(t) — p{,\) (filled curve) onto the 6^-axis. 
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We focus on the Mn2a-Mnl-Mn2b bond, involving 
only oxygens in the basal plane of the Mn2-pyramids. 
As illustrated in Fig.[lJ the corresponding spins in the C- 
phase of the real material have tii projections onto the 
a-direction. This is also the situation described by our 
model calculation. The LSDA+U charge density and net 
spin density within the plane of the Mn2a-Mnl-Mn2b 
bond are illustrated in Figs.|4](a) and (b). The densities 
in these figures correspond to the valence electrons within 
1.5 eV of the Fermi energy. As can clearly be observed 
in these figures, the tii spin configuration corresponds 
to a spatial distribution of the valence charge density, 
which breaks the center of inversion [circle in Fig.|4](a)]. 
In particular, the hybridization within a Mn2-0-Mnl 
bond depends on the relative Mn-spin orientation. This 
spin-dependent hybridization has a large effect on the 
oxygen states and yields a ferroelectric moment along 
6, as observed experimentally. The RSXS signal at the 
O i^-edge, which measures the magnetic polarization of 
oxygen, is directly related to this spin-dependent Mn2- 
0-Mnl hybridization and, hence, refiects Ve\. 

The above calculations do not take into account spin 
orbit interactions or structural distortions. Thus the FE 
polarization in this model results purely from the spin- 
dependent Mn2-0-Mnl hybridization and is due entirely 
to a redistribution of valence charge. The atomic posi- 
tions are not relaxed and therefore, by definition, there 
is no Vion- The important result therefore is that, within 
LSDA+U, a purely electronic contribution to the FE po- 
larization exists, which involves a spatial modulation of 
the 2p valence states (Fig.|4|. 

We also checked the infiuence of structural distortions 
within LSDA+U and found that they also affect the Mn- 
O covalency. However, in our calculations this effect is 
smaller than that due to the spin-induced valence charge 
redistribution, because the changes in atomic positions 
are very small. They are estimated to be of the order of 
0.01 A or less [71 [15] and, in fact, have not be resolved 
experimentally to date [26] . 

The LSDA+U calculations presented here agree very 
well with a more detailed DFT study of the related mate- 
rial HoMn2 05 [7 , which also revealed a covalency- driven 
Ve\ and showed that this contribution partially cancels 
the ionic V{on- The present experiments are confirmation 
of the theoretical prediction and indicate that Ve\ is a 
general and important feature of these compounds. 

We remark that the fact that T^eb ^is measured in our 
experiments, and V have the same temperature depen- 
dence (Fig. 2) indicates that either Ve\ is the dominant 
contribution to or that Ve\ and V^on have the same 
temperature dependence since V= 7^ion+ From our 
experiments we cannot a priori determine which of these 
is true. 

In summary, we have performed RSXS experiments 



on YMn205, revealing a spatial modulation of the 2p 
states, which is tightly coupled to the magnetic order, 
displaying the same periodicity and undergoing changes 
at the magnetic transition. Most importantly, the inte- 
grated intensity of the oxygen superstructure refiection is 
directly proportional to the macroscopic FE polarization. 
The present experiments verify that Ve\ does indeed ex- 
ist in a real material and imply that Mn-0 covalency and 
the resulting Ve\ play a major role for the ferroelectricity 
in YMnsOs. 
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